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The Schwinger mechanism, the production of charged particle-antiparticle pairs in a macroscopic 
external electric field, is derived for 2+1 dimensional theories. The rate of pair production per unit 
area for four species of massless fermions, with charge q, in a constant electric field E is given by 
7r~ 2 h~ 3 ^ 2 (T 1 / 2 (qE) 3 ^ 2 where c is the speed of light for the two-dimensional system. To the extent 
undoped graphene behaves like the quantum field-theoretic vacuum for massless fermions in 2+1 
dimensions, the Schwinger mechanism should be testable experimentally. A possible experimental 
configuration for this is proposed. Effects due to deviations from this idealized picture of graphene 
are briefly considered. It is argued that with present day samples of graphene, tests of the Schwinger 
formula may be possible. 



I. INTRODUCTION 

The Schwinger mechanism refers to the production 
of charged particle-antiparticle pairs out of the vacuum 
by a classical electric field which is homogenous over a 
large volume in space. More than half a century ago, 
Schwinger computed [l[ the rate associated with vacuum 
breakdown via pair production per unit volume for the 
case of a constant electric field. Over the years, the 
Schwinger mechanism has spawned a vast literature. In- 
voked to gain insights on topics as diverse as the string 
breaking rate in and on black hole physics @, 

this mechanism has become a textbook topic in quantum 
field theoryH. Topics such as back reaction || and finite 
size effects @] have been addressed. 

A key issue is the rate of pair production from a static 
electric field. Schwinger originally addressed this ques- 
tion by noting that the probability that a system with 
zero electrons or positrons remains in the fermionic vac- 
uum state decays exponentially in time[l|. Treating the 
electric field classically — i.e. the formal limit of q — > 0, 
E — > oo with qE fixed — Schwinger calculated the vacuum 
persistence probability, P vac (t), as a function of time: 

P vac (t) ee |(vac|C/(t)|vac)| 2 = exp(-wVt) (1) 
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where V is the spatial volume of the system and w is the 
rate of vacuum decay per unit volume. 

Schwingcr's interpretation[l[ of Eqs. ([T]) and © was 
straightforward: w was taken to be the local rate of pro- 
duction per unit volume of fermion-antifermion pairs by 
the electric field. This interpretation has been widely ac- 
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ccptcd in much of the literature on the Schwinger mech- 
anism. However, while the Schwinger formula of Eq. ([2|) 
is very well known, it has been argued that its inter- 
pretation as the pair production rate is not correct 12 1. 
Despite the very natural interpretation of w in Eq. ([2]) 
as the rate of production of pairs per unit volume, an 
explicit calculation gives the rate of pair production per 
unit volume, which we denote T, as 



r = 



(qE) 2 
An 3 h 2 c 



exp 



qEh 



(3) 



This rate does not agree with w: the entire rate is given 
by the first term in the series for w. For a recent, peda- 
gogical, discussion highlighting the theoretical distinction 
between these two rates, see Rcf. 
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Despite its theoretical significance, there has been no 
direct experimental signature of the Schwinger mecha- 
nism, of charged pair creation in electric fields. This 
is particularly unfortunate given the common confusion 
between the rate of pair creation T and w the rate of 
vacuum decay. It would be very useful to concoct an ex- 
perimental test to distinguish between the two directly. 
Moreover, apart from the distinction between w and T 
the derivation of the rate of pair production via the 
Schwinger mechanism raises a number of subtle issues as- 
sociated with the implementation of appropriate bound- 
ary conditions [3. 15, M. flp| : it is important to test whether 
these are handled correctly. Ultimately the most com- 
pelling test would be experimental. 

The reason that the Schwinger mechanism has never 
been tested experimentally is very easy to understand: 
the exponential factor in the pair production rate is very 
small for static macroscopic E fields realizable in the 
lab. It only becomes of order unity when E is large 
enough so that qE times the electron's Compton wave- 
length is greater than mc 2 : this requires an electric field 
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of order 10 16 V/cm; for E 

exp (-4 x 10 10 ). 

One might hope to test the Schwinger formula ex- 
perimentally in a condensed matter system which sim- 
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ulates light or masslcss, electrically charged, rclativis- 
tic fermions. In this context, the condensed matter 
system acts as an analog computer to test the under- 
lying result from relativistic field theory. Fortunately, 
it has been known for more than two decades that 
charged quasi-particle excitations in a potential with a 
two-dimensional hexagonal symmetry have a region of 
momenta over which their dispersion relation is lincar- 
(e — e ) 2 = c 2 (pl + p 2 )[ll[. This is precisely the dis- 
persion relation of a massless relativistic particle with 
energy measured relative to eo an d has c playing the role 
of c. Graphene (i.e., a single sheet of graphite) has such 
a symmetry. Moreover, in undoped graphene, the Fermi 
level is at eo. Thus, to the extent that a single parti- 
cle description holds in graphene, the quantum ground 
state of a filled Fermi sea is the precise analog of a filled 
Dirac sea — i.e. the vacuum of a two-dimensional non- 
interacting field theory for fermions. The recent devel- 
opment of techniques to produce samples of graphene 
and measure its properties [13, [l5[ has focused signifi- 
cant attention to its analogy with massless Dirac par- 
ticles: graphene has been proposed as a testing ground 
for the standard relativistic quantum mechanical effects 
of zwitterbewegung and Klcin/Landau-Zencr tunneling 
EH- This letter explores the possibility of using 
graphene to test experimentally the more subtle dynam- 
ics of the Schwinger mechanism. 



II. THE SCHWINGER PAIR CREATION RATE 
IN GRAPHENE 

One can adapt Schwinger's calculation [l| for smaller 
dimensions [tI M. [lOj ; in (2+1) dimensions the probability 
that the system has remained in the (fcrmionic) vacuum 
after time t, P vac (t), is 



P 2 +\t) = exp(-w 2+1 At) 
f(qE?' 2 



(4) 



with w 2+1 
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where A is the spatial area, c is the speed of light for the 
2-d system, £ is the Ricmann zcta function with £(3/2) 
2.612, and / is the number of species of fermion (i.e., four 
for graphene) . Similarly, the local rate of pair creation, 
r 2+1 , is given by the first term of this series fToL [Til. ll2T] . 
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These derivations hold for classical, constant, and exter- 
nally fixed electric fields. Consequently these relations 
are valid only when it is legitimate to neglect: i) macro- 
scopic back reaction to the applied field due to this same 
charged particle production rateQ; ii) the production of 
real photons as the charged particles accelerate in the 




FIG. 1: The dispersion relations are shifted on either side of 
the electric field. The image on the left is a cartoon picture of 
the in-state vacuum, where the Schwinger mechanism applies. 
After Tfm the level occupations change, and the system falls 
into a p-n type system, depicted on the right. 



electric field; and iii) interactions between the fermions 
mediated by the exchange of virtual photons. 

Alternative derivations [H [TTJ] to Schwinger's original 
calculation [l|, focus on the nature of a single particle level 
for the Dirac equation in the presence of an electric field, 
switched on in the distant past. These derivations have 
the twin virtues of clarifying both the role of boundary 
conditions, and the distinction between the rate associ- 
ated with vacuum breakdown, w, and the rate of pair 
creation, T. Here we will discuss the derivation done in 
the time independent gauge, Ei = — <9jj4o (ajj ) [3, [g]. In 
this gauge, we see that the electric field alters the sin- 
gle particle energy levels. The shifts in energy due to 
the potential have opposite signs on either side of the 
field region, allowing filled levels on one side to become 
degenerate with empty ones on the other. This yields 
an effective potential which filled levels in the elevated 
Dirac sea tunnel through, leaving a hole on one side and 
yielding a particle type state on the other. 

This tunneling problem is conceptually simple for an 
electric field of limited spatial extent Q. Consider an in- 
finite plane with an electric field independent of y, with 
magnitude E, oriented in the — x direction and confined 
to the region between —L/2 < x < L/2. A useful basis is 
the in-state wave functions which correspond to solutions 
of the Dirac equation with unit flux moving towards the 
region of the electric field from either the left or right. 
The states have amplitude T (times a flux-normalizing 
kinematic factor) to be found on the far side of the field 
region. T may be computed directly as a tunneling prob- 
lem in an effective energy-dependent Schrodingcr equa- 
tion derivable from the underlying Dirac equation. 

Suppose that the system is in the "in-state vacuum" 
(the left image in Fig. [I}. All of the the in-states below 
the Dirac sea on the left, ip^^,, are occupied for e < 
—qEL/2 and all of the states below the Dirac sea on the 
right, ip^p^, are occupied for e < qEL/2. Turning on the 
electric field in the distant past has merely shifted the 
energies of the states relative to the local vacuum, not 
their occupation number. The pair production rate for 
a particle with energy e and transverse momentum px is 
proportional to the transmission probability for a filled 
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in-state from the left —qEL/2 < e < qEL/2. In essence 
a filled level moving in the elevated Dirac sea towards 
the region of the electric field propagates through and 
emerges on the other side, with probability |T| 2 , where 
it appears as a particle. The rate of pair production per 
unit width per unit e per unit pr is computed analogously 
to the 3+1 dimensional caseQ. Integration over e and 
Pt gives the rate of pairs per unit width: 
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In the large L limit, the WKB result is increasingly valid, 



and \T\ [ 
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) Evaluating T within 
the WKB approximation for arbitrary L, and equating 
f t J^y with T2+1, one arrives at the 2+1 dimensional 
Schwinger formula of Eq. © for L 2 > Tic/qE. If L is 
not large one must compute T from the Dirac equation 
and numerically integrate over e and px- The ratio of 
the rate at finite L to the Schwinger rate depends only 
on the ratio L/ y/hc/qE. The numerically calculated rate 
smoothly converges to the Schwinger rate; see Fig. O 

It should be noted that while the preceding derivation 
is modeled on the derivation (in 3+1 dimensions) of Ref. 
Q there is one crucial distinction: Ref. @ computed T 
(the vaccum decay rate) which was implicitly assumed 
to be the pair production rate w. This difference only 
appears in the last line of Eq. ((7|): to compute W2+1 
rather than T 2 +i one replaces — log(l — \T(c,pt)\ 2 ) (used 
in Ref. i) with \T{e, PT )\ 2 . 



III. A PROPOSED EXPERIMENT 

The previous derivation depends on the system being 
in the in-state vacuum. Intuitively, after a certain tran- 
sient time during which the system equilibrates, the rate 
should be dominated by incoming levels from far away. 
Since the time scale for such transient behavior is finite 
at infinite L^, the natural time scale associated with 
transients does not depend on L. This transient time 
scale can be obtained via simple dimensional analysis: 

Ttrans = \/Tl/{cqE). 

The analysis crucially requires the fermionic "vacuum" 
to be placed in an external electric field which remains 
constant over time. To experimentally realize this with 
graphene, one must fix an external voltage by placing 
a region of the graphene sheet between two conducting 
plates, held at constant voltage. Pair production is driven 
by differences in level occupation on either side of the 
field. Particles and holes created in the field region are 
carried into the regions of the graphene sample on ci- 
ther side of the field regions which serve as reservoirs for 
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FIG. 2: The ratio of the pair production rate at finite L to 
the (infinite L) Schwinger formula rate for a 2+1 dimensional 
system. L is measured in units of Lo = \jTic/ (qE). 



particles and holes. After a finite time, the accessible 
states fill substantially and the system leaves the regime 
of validity of the Schwinger formula; the reservoirs de- 
velop an excess of holes or particles and the system more 
closely resembles p-n junctions, under extensive study in 
the context of Landau-Zener tunneling[l8l|; see Fig. [1] 
The natural timescale needed to deplete a substantial 
fraction of the reservoir and depart from the Schwinger 
regime is r m = ^qE/(hc 3 ) LL ICS . 

An experiment to test the Schwinger formula is thus 
conceptually straightforward: a sheet of graphene of 
width W, and length 2L ros + L is placed in an apparatus 
whose cross-section is given schematically in Fig. [3J L is 
the length the field region and L res is the length of the 
reservoirs to either side. An electric field of magnitude 
Vq/L is turned at t = 0; the Schwinger formula © should 
be accurate for Tt rans <C t <C T n n, bringing about a two- 
dimensional current density J from the Schwinger pairs. 
For r trans Ct< Tfin, the current density just beyond the 
field regions to good approximation is 



F>2+1 



Jsch = gr 



L = 



,3/2 



7T 2 H 3/2 C^L 1 / 2 



(8) 



As charge flows into the reservoirs, an external current 
/ = JschW^ must flow to the conductors to maintain them 
at fixed voltages of ±Vb/2. This current / can be mon- 
itored to determine J. Note that the graphene sheet 
is fully insulated electrically and is not part of a closed 
circuit. Thus, the current is necessarily transient and 
the experiment does not measure the usual conductance. 
This proposed method of testing Schwinger mechanism 
in graphene has a key feature in common with conduc- 
tance in graphene p-n junctions (GPNJ) in the ballistic 
regime: both systems depend on quantum tunnelling of 
massless Dirac particles jl8j. 

Effects due to non-zero temperature, non-zero lattice 
spacing, impurities, finite size effects and temporal tran- 
sients present in real systems can affect the results of the 
measurements. One expects to be in the regime of valid- 
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FIG. 3: Schematic depiction of the cross-sectional view of 
a possible experiment measuring the rate of production of 
Schwinger pairs. 



ity for the Schwinger formula if the parameters satisfy: 




T < {^U{qE a fL 2 /k 
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(9) 
(10) 

(11) 
(12) 



where Eq = Vq/L, a is the lattice spacing and l m f p is 
an effective mean-free path. Conditions ([9]) and (fT0|) for 
L relate to time transient effects and finite length is- 
sues associated with the WKB formalism, respectively; 
they are discussed above. The analogous condition for 
W follows from the requirement that the discrete mode 
sum in pt is approximated by the Gaussian integral in 
Eq. ([7]). Since the distance scale for the creation of 
pairs in the Schwinger mechanism is Lq = yjTic/{qE), 
the Schwinger formula only applies when the dynamics 
are well described by the simple Dirac description over 
that scale, requiring l m f p S> Lq. Condition (jll|l encodes 
the requirement of a linear dispersion relation; this fails 
for momenta comparable to a. Finally, Condition (fl"2"|) 
applies to thermal fluctuations, which imply a density 
of particles and holes in the reservoirs. These can ran- 
domly wander into, and get transported across, the field 
region, creating a current. The condition ensures that 
the Schwinger current dominates. 

The experimental configuration does not directly mea- 
sure conductivity since current flow is necessarily tran- 
sient. Still one might worry that the dynamics associated 
with the usual conductivity could mask the Schwinger 
effect. However, Eqs. © and © imply that Js c h = 



4rt rj_i 

h / L2ttJ 



L 2 qE 



. The term in parenthesis is of the 

scale expected from standard conductivity mechanisms: 
in graphene it is a ~ 4q 2 //i[l5j. The term in square 
brackets is of order unity; Condition (flT))) implies that the 



square root factor is large. Thus, the transient current 
density induced by the Schwinger mechanism dominates 
over what is expected from the usual conductivity. 

Whether Conditions ©-(H^J) can be satisfied in prac- 
tice depends critically on both the size of the graphene 
sample and its purity. In practice, with samples with 
sizes restricted to ~ 100 /xm and impurity concentrations 
reported in [U HH , the conditions appear to require ex- 
tremely fast measurements but do not seem to be beyond 
current technology. For the purpose of making estimates 
we will assume that the size of the sample is <~ 100 /im; 
for concrcteness we will take W = 100 /im, L = 1 /im 
and I/ res = 49 /im. We also take Vq = IV which ef- 
fectively satisfies Condition (|11|) [13]. With the values 
above, y/hc/qEo = 25nm, and Condition (fT0|) is satisfied 
by a factor of 40 for L and a factor of 4000 for W . Con- 
dition {T2J) is also well satisfied provided T < 1800°K. 

Condition (jTUJ) for l m f p depends on a quasi-particlc's 
energy, e. It has been argued with current samples the 
dominant contribution is from Coulomb impurities [l^]. 
For the purposes of estimating l m fp, this will be assumed 
to be correct. Coulomb scattering has an infinite cross- 
section and the mean-free path is not well defined. How- 
ever, an effective mean free path in the sense of the char- 
acteristic distance a quasiparticle travels before it is sub- 
stantially affected by the impurities can be estimated. 
It is the distance a quasi-particlc travels before enter- 
ing a region in which the Coulomb energy is comparable 
to the kinetic energy, yielding l m f p ~ Ke(?icani m p) _1 
where k the dielectric constant of the insulator and 
n-imp is the density of impurities. For Schwinger pairs, 
e is of order (but less than) qVo; Condition (fT0|) be- 
Ikl oo — w. m . as small 



comes 
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hcarii 

as ~ 2 x 10 11 cm -2 have been reported [ljS|]. Using this 
value, the left-hand side of the inequality is 350 k ran. 
The right-hand side is 25 nm. Thus, the condition ap- 
pears to be satisfied at least moderately well. Using the 
results in Fig. [2] as a guide in estimating errors sug- 
gests that at the least, a semi-quantitative test of the 
Schwinger mechanism should be possible. If a substrate 
with relatively large k proves viable, the condition may 
be well-satisfied. It appears possible that the mean-free 
paths are long enough to allow for a meaningful test of 
the Schwinger mechanism with currently available sam- 
ples. Precision tests will probably require cleaner samples 
which one hopes may become available in the future. 

With the parameters given above, 
1.9 x 10~ 9 s while T trans e 



TfiU 



qVoL j 



hL 

qV c 



2.6 x 10~ 14 s. It is easy to ensure t ^> T trans . The 
restriction t <C Tfm, however, requires taking data at 
a very high rate — considerably faster than one GHz. 
Fortunately, it is possible to take data at rates much 
faster than one GHz. If future sample sizes increase 
significantly, one could increase the size of L res and 
thereby Tfm and thus reduce the technical challenges 
associated with very rapid measurements. 
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IV. CONCLUSION 

To summarize, even with presently available samples, 
there is good reason to believe that the regime of valid- 
ity of the Schwingcr formula can be realized experimen- 
tally, at least at a semi-quantitative level. As larger and 
higher quality samples become available in the future, 
practical tests of the Schwinger formula with increasing 
accuracy ought to become possible. It is reasonable to 
expect that such experimental probes should become suf- 
ficient to test quantitatively the pair production of the 
Schwinger mechanism; it is is important that such mea- 



surements are accurate enough to distinguish between 
the predicted rate of pair of production and the rate of 
vacuum decay. 
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